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Abstract 

We present validation and results from a simulation of the Fermilab Tevatron including multiple 
beam dynamics effects. The essential features of the simulation include a fully 3D strong-strong 
beam-beam particle-in-cell Poisson solver, interactions among multiple bunches and both head-on 
and long-range beam-beam collisions, coupled linear optics and helical trajectory consistent with 
beam orbit measurements, chromaticity and resistive wall impedance. We validate individual phys- 
ical processes against measured data where possible, and analytic calculations elsewhere. Finally, 
we present simulations of the effects of increasing beam intensity with single and multiple bunches, 
and study the combined effect of long-range beam-beam interactions and transverse impedance. 
The results of the simulations were successfully used in Tevatron operations to support a change 
of chromaticity during the transition to collider mode optics, leading to a factor of two decrease in 
proton losses, and thus improved reliability of collider operations. 
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I. MOTIVATION 

The Fermilab Tevatron [l| is a p-p collider operating at a center-of-mass energy of 1.96 TeV 
and peak luminosity reaching 3.53 x 10'^^ cm^^ s^^. The colliding beams consist of 36 bunches 
moving in a common vacuum pipe. For high-energy physics operations, the beams collide 
head-on at two interation points (IPs) occupied by particle detectors. In the intervening 
arcs the beams are separated by means of electrostatic separators; long-range (also referred 
to as parasitic) collisions occur at 136 other locations. Effects arising from both head-on 
and long-range beam-beam interactions impose serious limitations on machine performance, 
hence constant efforts are being exerted to better understand the beam dynamics. Due to 
the extreme complexity of the problem a numerical simulation appears to be one of the most 
reliable ways to study the performance of the system. 

Studies of beam-beam interactions in the Tevatron Run II mainly concentrated on the 
incoherent effects, which were the major source of particle losses and emittance growth. This 
approach was justified by the fact that the available antiproton intensity was a factor of 10 
to 5 less than the proton intensity with approximately equal transverse emittances. Several 
simulation codes were developed and used for the optimization of the collider performance 

B.Q. 

With the commissioning of electron cooling in the Recycler, the number of antiprotons 
available to the collider substantionally increased. During the 2007 and 2008 runs the initial 
proton and antiproton intensities differed by only a factor of 3. Moreover, the electron 
cooling produces much smaller transverse emittance of the antiproton beam (~ 47rmmmrad 
95% normalized vs. ~ 207rmmmrad for protons), leading to the head-on beam-beam tune 
shifts of the two beams being essentially equal. The maximum attained total beam-beam 
parameter for protons and antiprotons is 0.028. 

Under these circumstances coherent beam-beam effects may become an issue. A num- 
ber of theoretical works exist predicting the loss of stability of coherent dipole oscillations 



when the ratio o: 
Landau damping 



Deam-beam parameters is greater than ~ 0.6 due to the suppression of 
4]. Also, the combined effect of the machine impedance and beam-beam 
interactions in extended length bunches couples longitudinal motion to transverse degrees 
of freedom and may produce a dipole or quadrupole mode instability ^. 

Understanding the interplay between all these effects requires a comprehensive simulation. 
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This paper presents a macroparticle simulation that includes the main features essential for 
studying the coherent motion of bunches in a collider: a self-consistent 3D Poisson solver 
for beam-beam force computation, multiple bunch tracking with the complete account of 
sequence and location of long-range and head-on collision points, and a machine model 
including our measurement based understanding of the coupled linear optics, chromaticity, 
and impedance. 

In Sections [UHV] we describe the simulation subcomponents and their validation against 
observed effects and analytic calculations. Section |Vl] shows results from simulation runs 
which present studies of increasing the beam intensity. Finally, in Section IVIII we study the 
coherent stability limits for the case of combined resistive wall impedance and long-range 
beam-beam interactions. 



II. BEAMBEAM3D CODE 

The Poisson solver in the BeamBeam3d code is described in Ref. 61]. Two beams are 
simulated with macroparticles generated with a random distribution in phase space. The 
accelerator ring is conceptually divided into arcs with potential interaction points at the ends 
of the arcs. The optics of each arc is modeled with a 6 x 6 linear map that transforms the 
phase space {x, x', y, y', z, 6} coordinates of each macroparticle from one end of the arc to the 
other. There is significant coupling between the horizontal and vertical transverse coordi- 
nates in the Tevatron. For our Tevatron simulations, the maps were calculated using coupled 



lattice functions 



a 



of beam element configuration to beam 



obtained by fitting a model 
position measurements. The longitudinal portion of the map produces synchrotron motion 
among the longitudinal coordinates with the frequency of the synchrotron tune. Chromatic- 
ity results in an additional momentum-dependent phase advance Sfix(y) = fioCx(y)Ap/p where 
Cx{y) is the normalized chromaticity for x (or y) and fiQ is the design phase advance for the 
arc. This is a generalization of the definition of chromaticity to apply to an arc, and reduces 
to the normalized chromaticity {Av/v)/{/S.p/p) when the arc encompasses the whole ring. 
The additional phase advance is applied to each particle in the decoupled coordinate basis 
so that symplecticity is preserved. 

The Tevatron includes electrostatic separators to generate a helical trajectory for the 
oppositely charged beams. The mean beam offset at the IP is included in the Poisson field 
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solver calculation. 

Different particle bunches are individually tracked through the accelerator. They interact 
with each other with the pattern and locations that they would have in the actual collider. 

The impedance model applies a momentum kick to the particles generated by the dipole 
component of resistive wall wakefields [9,]. Each beam bunch is divided longitudinally into 
slices containing approximately equal numbers of particles. As each bunch is transported 
through an arc, particles in slice i receive a transverse kick from the wake field induced by 
the dipole moment of the particles in forward slice j: 

Apl _ 2 / ATTeoc Njro<fj> L 

The length of the arc is L, Nj is the number of particles in slice j, tq is the classical 
electromagnetic radius of the beam particle e^/47reomoC^, % is the longitudinal distance 
between the particle in slice i that suffers the wakefield kick and slice j that induces the 
wake, fj is the mean transverse position of particles in slice j, 6 is the pipe radius, c is 
the speed of light, a is the conductivity of the beam pipe and are Lorentz factors of the 
beam. Quantities with units are specified in the MKSA system. 



III. SYNCHROBETATRON COMPARISONS 



We will assess the validity of the beam-beam calculation by comparing simulated synchro- 
betatron mode tunes with a measurement performed at the VEPP-2M 500 MeV e~^e~ collider 



and described in Ref. 



121 ]. These modes are an unambiguous marker of beam-beam inter- 



actions and provide a sensitive tool for evaluating calculational models. These modes arise 
in a colliding beam accelerator where the longitudinal bunch length and the transverse beta 
function are of comparable size. Particles at different z positions within a bunch are coupled 
through the electromagnetic interaction with the opposing beam leading to the development 
of coherent synchrobetatron modes. The tune shifts for different modes have a characteristic 
evolution with beam-beam parameter ^ = Nro/An'ye, in which is the number of particles, 
vq is the classical electromagnetic radius of the beam particle, and e is the unnormalized 
one-sigma beam emittance. 

There are two coherent transverse modes in the case of simple beam-beam collisions 
between equal intensity beams without synchrotron motion: the a mode where the two 
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beams oscillate with the same phase, and the vr mode where the two beams oscillate with 
opposite phases jlo|. Without synchrotron motion, the a mode mode has the same tune 
as unperturbed betatron motion while the vr mode frequency is offset by K^, where the 
parameter K is approximately equal to and greater than 1 and depends on the transverse 
shape of the beams pjj] . The presence of synchrotron motion introduces a more complicated 



spectrum of modes whose spectroscopy is outlined in Fig. 1 in Ref. [12|. 

We simulated the VEPP-2M collider using Courant- Snyder uncoupled maps. The hori- 
zontal emittance in the VEPP-2M beam is much larger than the vertical emittance. The 
bunch length (4 cm) is comparable to /?* = 6 cm so we expect to see synchrobetatron modes. 
In order to excite synchrobetatron modes, we set an initial y offset of one beam sigma 
approximately matching the experimental conditions. 

Longitudinal effects of the beam-beam interaction were simulated by dividing the bunch 
into six slices. At the interaction point, bunches drift through each other. Particles in 
overlapping slices are subjected to a transverse beam-beam kick calculated by solving the 
2D Poisson equation for the electric field with the charge density from particles in the 
overlapping beam slice. 

Simulation runs with a range of beam intensities corresponding to beam-beam parameters 
of up to 0.015 were performed, in effect mimicking the experimental procedure described 
in Ref. [l^. For each simulation run, mode peaks were extracted from the Fourier transform 
of the mean bunch vertical position. An example of the spectrum from such a run is shown 
in Fig. [1] with three mode peaks indicated. In Fig. [21 we plot the mode peaks from the 
BeamBeamSd simulation as a function of ^ as red diamonds overlaid on experimental data 



from Ref. [12] and a model using linearized coupled modes referred to as the matrix model 



described in that reference and Ref. 



13 



14| . As can be seen, there is good agreement 



between the observation and simulation giving us confidence in the beam-beam calculation. 

IV. IMPEDANCE TESTS 

Wakefields or, equivalently, impedance in an accelerator with a conducting vacuum pipe 
gives rise to well known instabilities. Our aim in this section is to demonstrate that the 
Wakefield model in BeamBeamSd quantitatively reproduced these theoretically and experi- 
mentally well understood phenomena. The strong head-tail instability examined by Chao 91] 
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FIG. 1: Simulated mode spectra in the VEPP-2M collider with ^ = .008 showing synchrobetatron 
modes. The line indicated by A is the base tune, B is the first synchrobetatron mode, C is the 
beam-beam tt mode. 

arises in extended length bunches in the presence of wakefields. For any particular acceler- 
ator optical and geometric parameters, there is an an intensity threshold above which the 
beam becomes unstable. 

The resistive wall impedance model applies an additional impulse kick in addition to the 
application of the map derived from beam optics. The tune spectrum is computed from the 
Fourier transform of the beam bunch positions sampled at the end of each arc. In order for 
the calculation to be a good approximation of the wakefield effect, the impedance kick should 
be much smaller than the x' or y' change due to regular beam transport so we divide the ring 
into multiple arcs, which brings up the question is how many is sufficient. The difference 
in calculated impedance tune shift for a 12 arc division of the ring or a 24 arc division 
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FIG. 2: The diamonds show simulated synchrobetatron modes as a function of beam-beam param- 
eter ^ (diamonds) and of observed modes (points). 

is only 2 x 10~^, which is less than 3% of the synchrotron tune (0.007 in this study), the 
relevant scale in these simulations. We perform the calculation with 12 arcs for calculational 
efficiency. 

In the absence of impedance, we would expect to see the tune spectrum peak at 20.574, 
the betatron tune of the lattice. With a pipe radius of 3 cm and a bunch length of 20 cm, 
resistive wall impedance produces the spectrum shown in Fig. [3] for a bunch of 4 x 10^^ 
protons at ISOGevf^]. In this simulation, the base tune z/^j is 20.574 and the synchrotron 
tune is 0.007. Three mode peaks are clearly evident corresponding to synchrobetatron modes 
with frequencies z/^ — Ug shifted up by the wakefield (point A), i/^ shifted down (point B), 
and z/^ + shifted upward (point C) as would be expected in Ref. jl5 |. 

In Fig. m we show the evolution of the two modes as a function of beam intensity. With 
the tune and beam environment parameters of this simulation, Chao's two particle model 
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FIG. 3: Simulated spectrum of a two slice bunch in the presence of wakefields and synchrotron 
motion showing three synchrobetatron modes A, B, and C induced by wakefields. 

predicts instability development at intensities of about 9 x 10^^ particles, which is close to 
where the upper and lower modes meet. We show two sets of curves for two slice and six 
slice Wakefield calculations. The difference between the two slice and six slice simulations 
is accounted for by the effective slice separation, z, that enters Eq. [H With two slices, 
the effective z is larger than than the six slice effective z, resulting in a smaller Wq. With 
the smaller wake strength, a larger number of protons is necessary to drive the two modes 
together as is seen in Fig. HI 

When the instability occurs, the maximum excursion of the bunch dipole moment grows 
exponentially as the beam executes turns through the accelerator. The growth rate can be 
determined by reading the slope of a graph of the absolute value of bunch mean position as a 
function of turn number plotted on a log scale. The growth rate per turn of dipole motion at 
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FIG. 4: Evolution of the base tune and lower synchrobetatron mode frequencies as a function of 
beam intensity showing the two modes approaching a common frequency due to impedance. The y 
scale is in units of the synchrotron tune. The simulations are shown for a two slice aond six slice 
Wakefield calculation. 

the threshold of strong head-tail instability has a parabolic dependence on beam intensity. 
The Wakefield calculation reproduces this feature, as shown in Fig. [51 The growth rate is 
slowly increasing up to the instability threshold at 5.42 x 10^^, after which it has the explicitly 
quadratic dependence on beam intensity (/) of growth rate = — 0.100 + 0.03041 — 0.00207/^. 

Chromaticity interacts with impedance to cause a different head-tail instability. We 
simulated a range of beam intensities and chromaticity values. The two particle model and 
the more general Vlasov equation calculation [ol indicate that the growth rate scales by the 
head-tail phase x = ^TrCvpz/ ct], where rj is the slip factor of the machine and z is roughly 
the bunch length. The head-tail phase gives the size of betatron phase variation due to 
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FIG. 5: The growth rate of dipole motion in the simulated accelerator with impedance as a function 
of beam intensity as the strong head-tail threshold is reached superimposed with a parabolic fit. 

chromatic effects over the length of the bunch. 

Some discussion of the meaning of the slip factor in the context of a simulation is neces- 
sary. In a real accelerator, the slip factor has an unambiguous meaning: rj = [ac — 1/7^)- 
The momentum compaction parameter ac is determined by the lattice and 7 is the Lorentz 
factor. We simulate longitudinal motion by applying maps to the particle coordinates z 
and 6 in discrete steps. The simulation parameters specifying longitudinal transport are the 
longitudinal beta function (3z and synchrotron tune Ug. Note that these parameters do not 
make reference to path length travelled by a particle. However, path length enters into the 
impedance calculation because wake forces are proportional to path length. In addition, 
analytic calculations of the effect of wake forces depend on the evolution of the longitudinal 
particle position which in turn depend explicitly on the slip factor. For our comparisons with 
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FIG. 6: The normalized growth rate of dipole motion in the simulated accelerator with impedance 
and chromaticity as a function of head-tail phase x three beam intensities demonstrating their 
linear relationship close to and the near- universal relationship for head-tail phase between — 1 
and 0. 

analytic results to be meaningful, we need to use a slip factor that is consistent with the lon- 
gitudinal maps and the path lengths that enter the wake force calculations. The relationship 
between the slip factor t] and the simulation parameters is (3z = 'nLr)/27ii>^, where Lq is the 
length of the accelerator and (3z = o'z/o's is the longitudinal beta function^, which may 
be derived by identifying corresponding terms in the solution to the differential equations 
of longitudinal motion and a one term linear map. 

When the growth rate is normalized by NroWo/27r(3'yi'j3, which includes the beam in- 
tensity and geometric factors, we expect a universal dependence of normalized growth rate 
versus head-tail phase that begins linearly with head-tail phase|3] and peaks around -1 23|. 
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Fig. [6] shows the simulated growth rate at three intensities with a range of chromaticites 
from —.001 to —0.5 to get head-tail phases in the to —1 range. The normalized curves are 
nearly identical and peak close to head-tail phase of unity. The deviation from a universal 
curve is again due to differences between the idealized model and detailed simulation. 



V. BUNCH-BY-BUNCH EMITTANCE GROWTH AT THE TEVATRON 

Understanding the effect of unwanted long-range collisions among multiple beam bunches 
in the design and operation of hadron colliders has received attention from other authors 
{2], 1^ which underscores the importance for this kind of simulation. A schematic of the 
fill pattern of proton and antiproton bunches in the Tevatron is shown in Fig. [71 There 
are three trains of twelve bunches for each species. A train occupies approximately 81.5° 
separated by a gap of about 38.5°. The bunch train and gap are replicated three times to 
fill the ring. Bunches collide head-on at the BO and DO interaction points but undergo long 
range (electromagnetic) beam-beam interactions at 136 other locations around the ring 24 1. 

Running the simulation with all 136 long-range IPs turns out to be very slow so we only 
calculated beam-beam forces at the two main IPs and and the long-range IPs immediately 
upstream and downstream of them. The transverse beta functions at the long-range collision 
locations are much larger than the bunch length, so the beam-beam calculation at those 
locations can be performed using only the 2D solver. 

One interesting consequence of the fill pattern and the helical trajectory is that any one 
of the 12 bunches in a train experiences collisions with the 36 bunches in the other beam 
at different locations around the ring, and in different transverse positions. This results 
in a different tune and emittance growth for each bunch of a train, but with the three- 
fold symmetry for the three trains. In the simulation, emittance growth arises from the 
effects of impedance acting on bunches that have been perturbed by beam-beam forces. 
The phenomenon of bunch dependent emittance growth is observed experimentally 201] . 



The beam-beam simulation with 36-on-36 bunches shows similar effects. We ran a simu- 
lation of 36 proton on 36 antiproton bunches for 50000 turns with the nominal helical orbit. 
The proton bunches had 8.8 x 10^^ particles (roughly four times the usual to enhance the 
effect) and the proton emittance was the typical 207rmmmrad. The antiproton bunch inten- 
sity and emittance were both half the corresponding proton bunch parameter. The initial 
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FIG. 7: Schematic of the position of proton and antiproton bunches in the Tevatron with 36 proton 
and 36 antiproton bunches. The diagram shows the positions at a time when the lead bunch of 
the trains are at the head-on colhsion location. Head-on collisions occur at location BO and DO. 
The green shading indicates the part of the ring where beam-beam collisions may occur in the 
simulations with six-on-six bunches. 
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FIG. 8: The simulated and measured emittance of each Tevatron proton bunch after running 
with 36 proton and 36 antiproton bunches. Curves (a) and (b) which show the emittance after 
50000 simulated turns are read with the left vertical axis. Curve (a) results from a simulation 
with the nominal beam spacing at the long-range IPs. Curve (b) results from a simulation with 
the hypothetical condition where the beam separation at the long-range IPs is 100 times normal, 
suppressing the effect of those long-range IPs. Curve (c) is the measured emittance of bunches 
after 15 minutes of a particular store {^5792) of bunches in the Tevatron, and is read with the 
right vertical axis. 



emittance for each proton bunch was the same so changes during the simulation reflect the 
beam-beam effect. 

Curve (a) in Figure [H] shows the emittance for each of the 36 proton bunches in a 36-on-36 
simulation after 50000 turns of simulation. The three-fold symmetry is evident. The end 
bunches of the train (bunch 1, 13, 25) are clearly different from the interior bunches. For 



14 



comparison, curve (c) shows the measured emittance taken during accelerator operations. 
The observed bunch emittance variation is similar to the simulation results. Another beam- 
beam simulation with the beam separation at the closest head-on IP expanded 100 times 
its nominal value resulted in curve (b) of Figure [S] showing a much reduced bunch-to-bunch 
variation. We conclude that the beam-beam effect at the long-range IPs is the origin of the 
bunch variation observed in the running machine and that our simulation of the beam helix 
is correct. 

VI. TEVATRON APPLICATIONS 

A. Single bunch features 

We looked at the tune spectrum with increasing intensity for equal intensity p and p 
beams containing one bunch each. As the intensity increases, the beam-beam parameter ^ 
increases. Fig. [9] shows the spectrum of the sum and difference of the two beam centroids 
for C, = 0.01,0.02,0.04, corresponding to beam bunches containing 2.2 x 10^^, 4.4 x 10^^ 
and 8.8 x 10^^ protons. The abscissa is shifted so the base tune is at and normalized in 
units of the beam-beam parameter at a beam intensity of 2.2 x 10^^. The coherent a and vr 
mode peaks are expected to be present in the spectra of the sum and difference signals of 
the two beam centroids. The coherent a modes are evident at 0, while the coherent tt modes 
should slightly greater than 1, 2, and 4 respectively. Increasing intensity also causes larger 
induced wake fields which broaden the mode peaks, especially the vr mode, as shown in 
Fig.l 

The 4D emittances at higher intensities show significant growth over 20000 turns as shown 
in Fig. [ini The kurtosis excess of the two beams remains slightly positive for the nominal 
intensity, but shows a slow increase at higher intensities indicating the the beam core is 
being concentrated as shown in Fig. [TTl Concentration of the bunch core while emittance is 
growing indicates the development of filamentation and halo. 

B. Simulation of bunch length, synchrotron motion and beam-beam interactions 

Synchrotron motion in extended length bunches modifies the effects of the beam-beam 
interaction by shifting and suppressing the coherent modes. The plots in Fig. [12] show 
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FIG. 9: Dipole mode spectra of the sum and difference offsets of two beam centroids at three beam 
intensities corresponding to beam-beam parameter values for each beam of 0.01, 0.02 and 0.04. 
The vertical scale is in arbitrary units. 

simulated spectra for sum and difference signals of the beam centroid offsets for one-on- 
one bunch collisions in a ring with Tevatron-like optics, with both short and long bunches, 
at three different synchrotron tunes. The sum signal will contain the a mode while the 
difference signal will contain the vr mode. In this Tevatron simulation, the beam strength is 
set so that the beam-beam parameter is 0.01, the base tune in the vertical plane is 0.576, and 
Py is approximately 30 cm. Subplots a and b of Fig. [12] show that with small synchrotron 
tune both the a and vr mode peaks are evident with short and long bunches. The a mode 
peak is at the proper place, with the vr mode peak shifted upwards by the expected amount, 
but with longer bunches (subplots c and d) the incoherent continuum is enhanced and the 
strength of the coherent peaks is reduced. When the synchrotron tune is the same as or 
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FIG. 10: The evolution of 4D emittances for beam-beam parameters of 0.01, 0.02, and 0.04 which 
correspond to intensities of (a) 2.2 x 10^\ (b) 4.4 x 10^\ (c) 8.8 x 10^\ and (d) 1.1 x 10^^ protons 
per bunch. 

larger than the beam-beam splitting (subplots e and /), short bunches still exhibit strong 
coherent modes, but with long bunches the coherent modes are significantly diluted. In the 
case of long bunches, the a mode has been shifted upwards to 0.580, and the vr mode is not 
clearly distinguishable from the continuum. At of 0.01 and 0.02, the synchrobetatron side 
bands are clearly evident. 

C. Multi-bunch mode studies 

When the Tevatron is running in its usual mode, each circulating beam contains 
36 bunches. Every bunch in one beam interacts with every bunch in the opposite beam, 
though only two interaction points are useful for high energy physics running. The other 
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turns 

FIG. 11: The evolution of (reduced) kurtosis of the particle distribution for intensities of (a) 2.2 x 
10^\ (b) 4.4 X 10", (c) 8.8 x lO", and (d) 1.1 x 10^^ protons per bunch. 

136 interaction points are unwanted and detrimental to beam lifetime and luminosity. The 
beam orbit is deflected in a helical shape by electrostatic separators to reduce the impact 
of these unwanted collisions, so the beams are transversely separated from each other in 
all but the two high-energy physics interaction points. Because of the helical orbit, the 
beam separation is different at each parasitic collision location. For instance, a bunch near 
the front of the bunch train will undergo more long-range close the the head-on interaction 
point, compared to a bunch near the rear of the bunch train. A particular bunch experiences 
collisions at specific interaction points with other bunches each of which has its own history 
of collisions. This causes bunch-to-bunch variation in disruption and emittance growth as 
will be demonstrated below. 

We will begin the validation and exploration of the multi-bunch implementation starting 
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FIG. 12: Simulated one-on-one bunch y plane a and vr mode tune spectra for short bunches (a, c, 
e) and long bunches (b, d, f), for three different synchtrotron tunes, with a Tevatron-like lattice. 

with runs of two-on-two bunches and six-on-six bunches before moving on to investigate the 
situation with the full Tevatron bunch fill of 36-on-36 bunches. Two-on-two bunches will 
demonstrate the bunches coupling amongst each other, but will not be enough to demon- 
strate the end bunch versus interior bunch behavior that characterizes the Tevatron. For 
that, we will look at six-on-six bunch runs. 

In these studies, we are only filling the ring with at most six bunches in a beam. Referring 
to Fig. [71 we see that only the head-on location at BO is within the green shaded region where 
beam-beam collisions may occur with six bunches in each beam. Because of the beam-beam 
collisions, each bunch is weakly coupled to every other bunch which gives rise to multi-bunch 
collective modes. 

We began the investigation of these effects with a simulation of beams with two bunches 
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each. The bunches are separated by 21 RF buckets as they are are in normal Tevatron 
operations. Colhsions occur at the head-on location and at parasitic locations 10.5 RF 
buckets distant on either side of the head-on location. To make any excited modes visible, 
we ran with 2.2 x 10^^ particles, which gives a single bunch beam-beam parameter of 0.01. 
There are four bunches in this problem. We label bunch 1 and 3 in beam 1 (proton) and 
bunch 2 and 4 in beam 2 (antiproton) with mean y positions of the bunches yi,...y4. 
By diagonalizing the covariance matrix of the turn-by-turn bunch centroid deviations, we 
determine four modes, shown in Fig. [131 Fig. [T3](a) shows the splitting of the a mode. The 
coefficients of the two modes indicate that this mode is primarily composed of the sum of 
corresponding beam bunches (1 with 2, 3 with 4) similar to the a mode in the one-on-one 
bunch case. The other two modes in Fig fT^ b) have the character and location in tune space 
of the 71 mode, from their coefficients and also their reduced strength compared to the a 
mode. 

With six on six bunches, features emerge that are clearly bunch position specific. 
Fig. fT4l(a) shows the turn-by-turn evolution of 4D emittance and (b) y kurtosis for each 
of the six proton bunches. It is striking that bunch 1, the first bunch in the sequence, has a 
lower emittance growth than all the other bunches. Emittance growth increases faster with 
increasing bunch number from bunches 2-5, but bunch 6 has a lower emittance growth than 
even bunch 4. The kurtosis of bunch 1 changes much less than that of any of the other 
bunches, but bunches 2-5 have a very similar evolution, while bunch 6 is markedly closer to 
bunch 1. One difference between the outside bunches (1 and 6) and the inside bunches (2-5) 
is that they have only one beam-beam interaction at the parasitic IP closest to the head-on 
collision, while the inside bunches have one collision before the head-on IP, and one after 
it. The two parasitic collision points closest to the head-on collision point have the smallest 
separation of any of the parasitics, so interaction there would be expected to disrupt the 
beam more than interactions at other parasitic locations. 

To test this hypothosis, we did two additional runs. In the first, the beam separation at 
the parasitic IP immediately downstream of the head-on IP was artificially increased in the 
simulation so as to have essentially no effect. The effect of this is that the first proton bunch 
will not have any beam-beam collisions at an IP close to the head-on IP, while all the other 
bunches will have one collision at a near-head-on IP. The corresponding plots of emittance 
and kurtosis are shown in Fig. [151 The kurtosis data shows that bunches 2-5 which all suffer 
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FIG. 13: Mode tune spectrum for a two on two bunch run at 2.2 X IQii particles/bunch = 0.01). 
Figure (a) shows the two modes that are most hke a modes, a mode 1 is 0.53yi + 0.53y2 + 0.59y3 — 
0.317/4, a mode 2 is 0.39yi + 0.49y2 — 0.46y3 — 0.63y4. Figure (b) shows the two 7r-hke modes, tt 
mode 1 is 0.74yi - O.662/2 - O.OSys, vr mode 2 is 0.12yi + 0.20y2 - O.662/3 + 0.31^4- The absolute 
scale of the Fourier power is arbitrary, but the relative scales of plots (a) and (b) are the same. 



one beam-beam collision at a close parasitic IP are all together while bunch 1 which does 
not have a close IP collision is separated from the others. 

Emittance and kurtosis growth in simulations where the beam separation at the closest 
upstream and downstream parasitic IPs was increased is shown in Fig. [161 In this configu- 
ration no bunch suffers a strong beam-beam collision at a parasitic IP close to the head-on 
location so the kurtosis of all the bunches evolves similarly. 
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FIG. 14: A six-on-six bunch Tevatron run with 8.8 x 10^^ particles/bunch: (a) The turn-by-turn 
evolution of 4D emittance of each of the six bunches, (b) The turn-by-turn evolution of y kurtosis 
of the six bunches. 
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FIG. 15: In a six-on-six bunch Tevatron run with 8.8 x 10^^ particles/bunch, with the beam 
spacing at the first parasitic IP downstream of the head-on location artificially increased: (a) The 
4D emittance of each of the six bunches as a function of turn, (b) the y kurtosis of the six bunches 
as a function of turn. 

VII. LOWER CHROMATICITY THRESHOLD 

During the Tevatron operation in 2009 the hmit for increasing the initial luminosity was 
determined by particle losses in the so-called squeeze phase [211]. At this stage the beams 
are separated in the main interaction points (not colliding head-on), and the machine optics 
is gradually changed to decrease the beta-function at these locations from 1.5 m to 0.28 m. 
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FIG. 16: In a six-on-six bunch Tevatron run with 8.8 x 10^^ particles/bunch, with the both nearest 
upstream and downstream parasitic IP artificially widened: (a) The 4D emittance of each of the 
six bunches as a function of turn, (b) the y kurtosis of the six bunches as a function of turn. 

With proton bunch intensities currently approaching 3.2 x 10^^ particles, the chromaticity 
of the Tevatron has to be managed carefully to avoid the development of a head-tail insta- 
bility. It was determined experimentally that after the head-on collisions are initiated, the 
Landau damping introduced by beam-beam interaction is strong enough to maintain beam 
stability at chromaticity of +2 units (in Tevatron operations, chromaticity is Au/{Ap/p).) 
At the earlier stages of the collider cycle, when beam-beam effects are limited to long-range 
interactions the chromaticity was kept as high as 15 units since the concern was that the 
Landau damping is insufficient to suppress the instability. At the same time, high chro- 
maticity causes particle losses which are often large enough to quench the superconducting 
magnets, and hence it is desireable to keep it at a reasonable minimum. 

Our multi-physics simulation was used to determine the safe lower limit for chromaticity. 
The simulations were performed with starting beam parameters listed in Table [H With 
chromaticity set to -2 units, and no beam-beam effect, the beams are clearly unstable as seen 
in Fig.[T71 With beams separated, turning on the beam-beam effect prevents rapid oscillation 
growth during the simulation as shown in Fig. [181 The bursts of increased amplitude is 
sometimes indicative of the onset of instability, but it is not obvious within the limited 
duration of this run. The RMS size of the beam also does not exhibit any obvious unstable 
tendencies as shown in Fig. [111 

Based on these findings the chromaticity in the squeeze was lowered by a factor of two. 
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TABLE I: Beam parameters for Tevatron simulation 



Parameter 



value 



beam energy 

p particles/bunch 

p particles/bunch 

p tune {i^x,'^y) 

p (normalized) emittance 

p tune {i^xji^y) 

p (normalized) emittance 

synchrotron tune Vs 

slip factor 

bunch length (rms) 

Sp/p momentum spread 

effective pipe radius 



980 GeV 
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3.0 X 10 
0.9 X 10^^ 
(20.585,20.587) 
20-7r mm mrad 
(20.577,20.570) 
Gtt mm mrad 
0.0007 
0.002483 
43 cm 
1.2 X 10-^ 
3 cm 



and presently is kept at 8-9 units. This resulted in a significant decrease of the observed 



particle loss rates (see, e.g.. Fig. 5 in |2l[|). 



VIII. SUMMARY 

The key features of the developed simulation include fully three-dimentional strong-strong 
multi-bunch beam-beam interactions with multiple interaction points, transverse resistive 
wall impedance, and chromaticity. The beam-beam interaction model has been shown to 
reproduce the location and evolution of synchrobetatron modes characteristic of the 3D 
strong-strong beam-beam interaction observed in experimental data from the VEPP-2M 
collider. The impedance calculation with macroparticles excites both the strong and weak 
head-tail instabilities with thresholds and growth rates that are consistent with expectations 
from a simple two-particle model and Vlasov calculation. Simulation of the interplay be- 
tween the helical beam-orbit, long range beam-beam interactions and the collision pattern 
qualitatively matches observed patterns of emittance growth. 

The new program is a valuable tool for evaluation of the interplay between the beam- 
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FIG. 17: The x dipole moment in a simulation with C = —2 and no beam-beam effect showing the 
development of instability. 

beam effects and transverse collective instabilities. Simulations have been successfully used 
to support the change of chromaticity at the Tevatron, demonstrating that even the reduced 
beam-beam effect from long-range collisions may provide enough Landau damping to prevent 
the development of head-tail instability. These results were used in Tevatron operations to 
support a change of chromaticity during the transition to collider mode optics, leading to a 
factor of two decrease in proton losses, and thus improved reliability of collider operations. 
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FIG. 18: The x dipole moment of a representative bunch in a 36-on-36 simulation with C = —2 
with beam-beam effects and beams separated showing no obvious instabihty within the hmits of 
the simulation. 
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FIG. 19: The x RMS moment of a representative bunch in a 36-on-36 simulation with C = —2 
with beam-beam effects and beams separated showing no obvious instabihty within the hmits of 
the simulation. 
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order to drive the strong headtail instability for comparison with the analytical model. 

The simulated machine is above transition (77 is positive.) The head-tail instability develops 

when chromaticity is negative, thus the head-tail phase is negative. 

With three-fold symmetry of bunch trains, train-on-train collisions occur at six locations 
around the ring. The collision of two trains of 12 bunches each results in bunch-bunch col- 
lisions at 23 locations which when multiplied by six results in 138 collision points. It is a 
straightforward computer exercise to enumerate these locations. Two of these locations are 
distinguished as head-on while the remainder are parasitic. [2C] 
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